Adiabatic passage for quantum gates in mesoscopic atomic ensembles 
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We present schemes for geometric phase compensation in adiabatic passage which can be used 
for implementation of quantum logic gates with atomic ensembles consisting of an arbitrary number 
of strongly interacting atoms. Protocols using double sequences of stimulated Raman adiabatic 
passage (STIRAP) or adiabatic rapid passage (ARP) pulses are analyzed. Switching the sign of the 
detuning between two STIRAP sequences, or inverting the phase between two ARP pulses, provides 
state transfer with well defined amplitude and phase independent of atom number in the Rydberg 
blockade regime. Using these pulse sequences we present protocols for universal single-qubit and 
two-qubit operations in atomic ensembles containing an unknown number of atoms. 

PACS numbers: 32.80.Ee, 03.67.Lx, 34.10.+X, 32.70.Jz , 32.80.Rm 
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Quantum information can be stored in collective states 
of ensembles of strongly interacting atoms [l| . This idea 
can be extended to encoding an entire register of qubits 
in ensembles of atoms with multiple ground states 
which opens up the possibility of large quantum registers 
in a single atomic ensemble 13|, or of coupling arrays of 
small ensembles in a scalable atom chip based architec- 
ture Quantum information based on ensembles can 
be realized more generally in any ensemble of strongly 
coupled spins 0. Our proposal for implementing high 
fidelity quantum gates in ensembles is thus of interest for 
several different implementations of quantum computing. 

The use of ensembles for qubit encoding is of interest 
for several reasons. The enhanced coupling to the ra- 
diation field by a factor of \/N, with the number of 
atoms or spins, is useful for coupling matter qubits to sin- 
gle photons [6|. Combining photon coupling with local 
quantum gates in ensembles enables architectures with 
improved fidelity for quantum networking [tI . The use of 
ensemble qubits is also attractive for deterministic load- 
ing of registers of single atom qubits |8| and for realizing 
gates that act on multiple particles [l3| • AH of these 
capabilities rely on high fidelity quantum gate operations 
between collectively encoded qubits. However, due to the 
increase of the Rabi frequency of oscillations between dif- 
ferent collective states proportional to a//V, ~ ^/Nfl, 
with r2 the one atom Rabi frequency, it is difficult to per- 
form gates with well defined rotation angles in the situa- 



tion where N is unknown jllL Il2| . Although there is re- 



cent progress in nondestructive measurement of N with 
high accuracy [l3j it remains an outstanding challenge 
to impleixient high fidelity quantum logic gates without 
precise knowledge of N, particularly in the case of col- 
lectively encoded registers Q where the effective value of 
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FIG. 1. (Color online), (a) Energy levels for two-photon STI- 
RAP and single-photon ARP excitation; (b) Time sequence 
of STIRAP laser pulses; (c) Time sequence for ARP laser ex- 
citation; (d) Time dependence of the single-atom excitation 
probability. 



A'^ depends on the unknown quantum state encountered 
during a computation. 

Adiabatic passage techniques have been widely used 
for deterministic population transfer in atomic and 
molecular systems [IJ, [15 1. These techniques have 
been studied for quantum state control Il6| . qubit ro- 
tations [l3], creation of entangled states lisl . and for 
deterministic excitation of Rydberg atoms (l9l [20| . Al- 



though STIRAP or ARP methods provide pulse areas 
with strongly suppressed sensitivity to the Rabi fre- 
quency JIat, and therefore suppressed sensitivity to TV, 
the phase of the final state is in general still strongly 
dependent on N . Increase of the average atom number 
NS reduces the relative fluctuations by a factor of \fN , 
which is still not enough for scalable quantum computa- 
tion with moderate size ensembles. Indeed gate errors at 
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FIG. 2. (Color online), (a) Poissonian statistics of atom 
loading into an optical trap for N = 5; (b) Comparison of 
the fidelity of single-atom excitation by a tt laser pulse with 
the area optimized for N — 5 atoms {t — -jr/y/EQ), with 
a STIRAP sequence, and an ARP pulse. The single STI- 
RAP sequence used D.j{t) = ^^6"''+*^^''/^^^ for j = 1,2 with 
Q.i/2tv = 30 MHz, Q.2/2tv = 40 MHz, ti = 3.5 ixs, t2 = 5.5 ^s, 
r = 1 /is, and S/2-k = 200 MHz. The single ARP pulse used 

no{t) = fioe"*''''^^^ with no/2-K = 2 MHz, r = lus, and 
linear chirp a/27r = (l/2TT){dS(t) /dt) = 1 MHz/^is [20]. 

the level of 0.001 would require N ^ 4000 ( accounting for 
fluctuations of ±2.5'\//V) , and achieving full blockade for 
such a large ensemble remains an outstanding challenge. 
In this Letter we propose double adiabatic sequences us- 
ing either STIRAP or ARP excitation which remove the 
phase sensitivity, and can be used to implement gates on 
collectively encoded qubits without precise knowledge of 
N. 

Our approach is shown in Fig. [TJ A sequence of two 
STIRAP pulses is produced with fields having Rabi fre- 
quencies ill, ^^2, and detuning 6 from the intermediate 
state. A Rydberg state \r) is excited by a counterintu- 
itive sequence of laser pulses in the regime of Rydberg 
blockade, where excitation of more than one Rydberg 
atom is suppressed due to strong dipole-dipole interac- 
tion between the atoms [l[ . The second reverse STIRAP 
sequence, as shown in Fig. [ija), returns the Rydberg 
atom back to the ground state. The time dependence of 
the probability of single-atom Rydberg excitation in an 
ensemble of two interacting atoms is shown in Fig. 1(d) 
for r^i = = 27r X 15 MHz and two different values 
of detuning. As we have shown in our recent work [20j . 
large detuning from the intermediate excited state |e) 
is required for deterministic single-atom excitation after 
the first STIRAP sequence, which is represented by the 
solid curve in Fig. [TJd). In the case of zero detuning, 
a single STIRAP sequence does not leave any popula- 
tion in the Rydberg state [broken curve in Fig. [ijd)], in 
agreement with the calculations of Ref. [l^l ■ Determin- 
istic single-atom excitation can also be achieved using 
linearly chirped ARP pulses, as shown in Figs. [IJc),(d). 

The probability of loading N atoms in a small optical 
or magnetic trap is described, in general, by Poissonian 
statistics. Figure 2(a) shows the Poissonian distribution 
for iV = 5. For iV = 5 the probability to load zero 
atoms is 0.0067, which is small enough to create a large 
quantum register with a small number of defects [2l|. 
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FIG. 3. (Color online). Calculated time dependence of 
the phase of the collective ground state amplitude for A*' = 
1, 2, 7 atoms (top to bottom). Double STIRAP sequence with 
S/2TT = 200 MHz (a), with S/2tv = 200 MHz x sign (t) (b), 
and for a double ARP pulse sequence with phase inversion 
(c). All other parameters are as in Fig. [2] 



Figure 2(b) shows a comparison of the fidelity of single- 
atom excitation for a single-photon tt rotation with the 
area optimized for N = 5 atoms compared to STIRAP or 
ARP pulses. We see that the adiabatic pulses reduce the 
population error by up to several orders of magnitude for 
a wide range of N. Calculations were performed using the 
Schrodinger equation, neglecting spontaneous emission, 
and assuming perfect blockade so only states with at most 
a single Rydberg excitation were included. 

At the end of a double STIRAP sequence the popula- 
tion is returned back to the collective ground state |000...) 
of the atomic ensemble, but a geometric phase is accumu- 
lated. This phase shift of the ground state is dependent 
on the Rabi frequency and leads to gate errors. The 
phase of the atomic wavefunction can be compensated 
by switching the sign of the detuning between two STI- 
RAP sequences, or by switching the phase between two 
ARP pulses, as shown in Fig. H For a double STIRAP 
sequence with the same detuning throughout the accu- 
mulated phase depends on [Fig. Efa)], while the phase 
change is zero, independent of N, when we switch the 
sign of detuning 6 between the two STIRAP sequences 
[Fig. El^b)]. A similar phase cancellation occurs for tt 
phase shifted ARP pulses [Fig. Efc)] Although the 

proposed double-pulse sequences are almost insensitive 
to variations of the Rabi frequency, they are not immune 
to errors caused by fiuctuations of the Rabi frequencies 
between the first and second pulses. The dependence of 
the population and phase errors on parameters of the 
laser pulses are shown in Fig. The population transfer 
error in the ensemble of A'^ = 5 atoms can be kept below 
10-3 for STIRAP [Fig.UJa)] and 10"* for an ARP pulse 
[Fig. m^b)] for a wide range of Rabi frequencies. The 
dependence of the phase error on the ratio of Rabi fre- 
quencies Q^f^ /fl']^^ between pulses [sec Fig.[IJb)] is shown 
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FIG. 4. (Color online). Dependence of the population error 
on the Rabi frequency with N = 5 atoms for STIRAP (a) 
or ARP pulses (b). Dependence of the phase error on Rabi 
frequency changes between pulses for STIRAP (c) or ARP 
pulses (d). 



in Fig. SJc) for N = 1..5 atoms. The ARP excitation in 
Fig. m^d) demonstrates less sensitivity to fluctuations of 
the Rabi frequency. For either approach the double pulse 
amplitudes must be well matched for low phase errors. 
Due to the fact that the pulses are separated in time by 
only a few us amplitude matching at the level of 10^^ is 



feasible [23 1. 



The phase compensated double STIRAP or ARP se- 
quences can be used to implement a universal set of 
quantum gates as we now describe. Consider atoms 
with levels |0), |1), |e)|?-) as shown in Fig. 1. A qubit 
can be encoded in an N atom ensemble with the logical 
states |0) = I000...000), |1)' = ^ Ej^i |000...1j...OOO). 
Levels |0),|1) are atomic hyperfine ground states with 
coupling between these states and implementation of 
quantum gates mediated by the singly excited Rydbcrg 
state |f)' = ^E^i |000-^i-000)- Rydberg block- 
ade only allows single excitation of |r) so the states |0) 
and |f)' experience a collectively enhanced coupling rate 
riN ~ \/Nfl. States |f)' and are coupled at the sin- 
gle atom rate il. State is produced by the sequential 
application of tt pulses |0) — )• |f)' and |f)' — )■ 

Pulse areas independent of N on the |0) O |r) tran- 
sition can be implemented with STIRAP or ARP as de- 
scribed above. We will define the logical basis states 
and the auxiliary Rydberg state as |0) = |000...000), 
|i) = e*'^"|l)', and |f) = e'>^"|f)'. Here xn is the phase 
produced by a single N atom STIRAP pulse with posi- 
tive detuning. We assume that we do not know the value 
of A^, which may vary from qubit to qubit, and therefore 
Xn is also unknown, but has a definite value for fixed N. 

For universal quantum computing we need a complete 
set of gates which can be comprised of single qubit rota- 
tions and a two-qubit entangling gate. Arbitrary single 




FIG. 5. (color online) (a) Single qubit gate for a mesoscopic 
qubit with A'" atoms. Pulses 1 — 5 act between the qubit 
states |0), 1 1) and the Rydberg states |ro), |ri). Pulses 1, 2, 4, 5 
are optical transitions and pulse 3 is a microwave frequency 
transition between Rydberg states, (b) CNOT gate between 
mesoscopic qubits with Nc atoms in the control qubit and Nt 
atoms in the target qubit. 



qubit rotations in the basis |0), |1) can be performed with 
high fidelity, without precise knowledge of A^, by access- 
ing several Rydberg levels as shown in Fig.[5Ia). Starting 
with a qubit state|'0) = a|0) -t-&|l) we perform a sequence 
of pulses 1-5, giving the sequence of states 
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The final state = a'|0) + b'\l) is arbitrary and is 



selected by the rotation R{d,(j)), in step 3: 



-b' 



R( 



Depending on the choice of implementa- 



tion, to be discussed below, this may be a one- or two- 
photon microwave pulse. Provided states |ro),|ri) are 
strongly interacting, and limit the entire ensemble to a 
single Rydberg excitation in either state, the indicated 
sequence is obtained. In the regime of ils large compared 
to the Rydberg excitation rates the time spent populat- 
ing a Rydberg level corresponds to 47r of Rydberg pulse 
area. This is the same as for a single atom Cz gate, 
and we therefore expect the limit on gate infidelity to 
be ^ 0.002 for small ensembles. For moderate size 
ensembles the gate error is expected to increase linearly 
with N [13]. The five pulse sequence we describe here 
is more complicated than the three pulses needed for an 
arbitrary single qubit gate in the approach of Ref. [l[. 
The reason for this added complexity is that the special 
phase preserving property of the double STIRAP or ARP 
sequences requires that all population is initially in one 
of the states connected by the pulses. The sequence of 
pulses in Fig. ^a) ensures that this condition is always 
satisfied. 

A CNOT gate can be implemented in several ways. 
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(2) 



The standard construction [25| of H(t) - Cz - H(t) 
can be used where the Hadamard gates are performed 
as in Fig. [5l^a), and the Cz operation is im plem ented 
in the same way as for single atom qubits [26| using 
7r|i)_|p)(c) 27r|j)_|p^ (t) 7r|j)_|p) (c), where c(t) stand for 
control(target) qulaits. The CNOT gate therefore re- 
quires a total pulse area of 12tt Rydberg pulses. We 
can reduce this to Ttt of Rydberg pulses as shown in 
Fig. [5l[b) which implements an approach analogous to 
the amplitude-swap gate demonstrated for single atom 
qubits in . Starting with an arbitrary two-qubit state 
= a|00) + 6|01) + c|10) -I- we generate the se- 

quence of states 

IV'i) = a|foO) + 6|foi) + c|10) + d\n) 
IV'2) = a|foO) + 6|rol) + c|10) + id\lfi) 
IV-s) = al^'oO) + b\fol) + c\lfo) + id\lfi) 

= m|fiO) -t- ib\fil) + ic\lfi) - d\lfa) 
1^5) = ia\fiQ) + ib\fil) + ic\lfi) - d|10) 
|^6> = ia\riQ) + ib\ril) - c\ll) - d\10) 
|^7> = ia\00) + ib\Ol) - c\ll) - d\10). 
The gate matrix is therefore 



C^CNOT 



which can be converted into a standard CNOT gate with 
a single qubit rotation. All pulses except 4 are optical 
and are localized to either the control or target qubit. 
Pulse 4 is a microwave field and drives a tt rotation on 
both qubits. As for the single qubit gate the requirement 
for correct operation is that the interactions |ro) o |ro), 
^ \ro) o \ri) all lead to full blockade of the 
ensembles. Since the frequency of pulse 4, which is de- 
termined by the energy separation of states \ro), |ri), can 
be chosen to be very different from the qubit frequency 
given by the energy separation of states |0), |1) the ap- 
plication of microwave pulses will not lead to crosstalk in 
an array of ensemble qubits. 

We now discuss the feasibility of implementation of 
these gate protocols in a Rydberg blockaded ensemble. 
To get an isotropic interaction suitable for ensemble 
blockade we use ns states The long range inter- 

action strength can be parameterized with a Cq coeffi- 
cient as V{n,n') = Cg"'" / with R the atomic sep- 
aration. For Cs ns states the optimum gate fidelity is 
obtained for 80s [2J|, and the interaction strengths for 
|ro) = |80si/2,m = 1/2), In) = |81si/2,m = 1/2) are 
Cr-''^ = 5.1, C^-''^ ^ 3.7, in units of 
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10® MHz ^m® . Rydberg ns states can be accessed start- 
ing from a ground s state using two-photon STIRAP 
pulses. These ns states can also be excited with two- 
photon ARP pulses where one photon is fixed frequency 



and one is chirped. Alternatively single photon ARP 
pulses can be used to access np states. Although the 
interaction of np states is not isotropic it can be made 
isotropic in lower dimensional 1- or 2-D lattices by ori- 
enting the quantization axis perpendicular to the lattice 
symmetry plane. For Cs the optimal state is 112p3/2 24| . 
and the interaction strengths for |ro) = |112p3/2i'^ = 
3/2), \ri) = |113p3/2,m = 3/2) at 90 deg. to the quan- 



tization axis are Cg^^^'^^^' = 250. 



a 



(112,113) 



820. 



(113,113) 



270., in units of 10*^ MHz /im*^. Wc see that 



for both ns and np states a strong interaction is obtained 
for all involved Rydberg states as desired. The pulse 
connecting |ro),|ri) can be implemented as a 2-photon 
electric dipole transition at microwave frequencies via 
a neighboring opposite parity state. The large transi- 
tion dipoles of Rydberg states scaling as n^eog (e is the 
electronic charge, ag is the Bohr radius) render fast mi- 
crowave pulses straightforward to implement. At n = 80, 
a detuning of 1 GHz from the intermediate state, and 
a very modest 1 /iW/cm^ microwave power level, gives 
~ 25 MHz two-photon Rabi frequency. 

The requirements for gates that are insensitive to the 
value of N are within reach of current experimental capa- 
bilities. Optimization of the Rabi frequency and a rigor- 
ous determination of the gate fidelity could be found fol- 
lowing the methods of Ref. . These collective gates re- 
quire more Rydberg pulses than single atom gates which 
will result in different choices for the Rabi frequencies 
and optimal quantum numbers n compared to the values 
used for illustration above which are optimal for single 
atom qubits. The influence of a finite strength Rydberg 
interaction must also be accounted for. Such a calcula- 
tion is beyond the scope of this Letter and will be pre- 
sented elsewhere. The collective gates do require that N 
does not change during the operation of the gate. In any 
realistic implementation with trapped atoms gate times 
are significantly shorter than atom loss times, so this is 
not an issue. On longer time scales atom loss leads to 
qubit errors at a physical level, as is also true for a single 
atom encoding. Ideas for correcting physical errors due 
to atom loss in ensembles have been presented in Ref. 

In summary wc have analyzed double STIRAP and 
ARP sequences with phase compensation for quantum 
gates in collectively encoded ensembles. Our analysis 
shows that high fidelity universal gates can be achieved 
using available experimental resources. We have made 
explicit calculations for N = 1..10 showing the perfor- 
mance of these pulse schemes. For larger ensembles, with 
smaller fractional variation in N, the performance will 
be even better. We anticipate that these ideas will con- 
tribute to realization of quantum logic using collectively 
encoded qubits and registers. 
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